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Abstract. 

The random motion of a Brownian particle confined in some finite domain is con- 
sidered. Quite generally, the relevant statistical properties involve infinite series, whose 
coefficients are related to the eigenvalues of the diffusion operator. Unfortunately, the 
C^ , latter depend on space dimensionality and on the particular shape of the domain, and 

an analytical expression is in most circumstances not available. In this article, it is 

, ■ shown that the series may in some circumstances sum up exactly. Explicit calculations 

^ , are performed for 2D diffusion restricted to a circular domain and 3D diffusion inside 

O ' a sphere. In both cases, the short-time behaviour of the mean square displacement is 

I ^i . obtained. 

^^ ■ Keywords. Brownian motion. Colloids. 
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1. Introduction 

Although the theory of Brownian motion has been one of the cornerstones of 
modern physics for more than one century |jy, it still raises numerous fundamental 
questions [2j. Applications are found in various disciplines, especially in biophysics where 
technical progress allows nowadays for the detection of individual nanoparticles in living 
systems j3]. In many instances, the particles' motion is strictly restricted to bounded 
domains. One can quote corralled motion of receptors on cell membrane jH El Ej, 
protein diffusion inside the cell nucleus j7], or transporter diffusion through nuclear 
pores [Hlin]- In a typical experiment, trajectories of tracer particles are recorded and 
analysed in terms of position correlations or mean square displacement. From those 
datas, information on the size of the diffusing objects or on their interactions with the 
surrounding medium can then be extracted |l()j . 

In this article, we consider the motion of a Brownian particle confined in a domain 
V of given size and shape. The statistical properties of the particle are described by 
the Green's function G (r, r', t), which represents the probability density of finding the 
particle at point r and time t, together with the intial condition G (r, r', 0) = 5 (r — r'). 
The Green's function satisfies the Fokker-Planck equation 

f)G 

^-DW'G = 5ir-r')5{t) (1) 

with D the diffusion coefficient. Confinement is enforced through the reflecting boundary 
conditions on the frontier dT> of the domain p. It should be noticed that, because of 
the finite size of the system, translational invariance is broken and the Green's function 
depends on r and r' separately. Once the propagator is known, all the statistical 
properties can, in principle, be evaluated. For instance, the mean square displacement 
(MSD) is obtained from 

{6r\t)) = {{r{t) - r{0)f) 

= f d'^r f dV G (r, r', t) Po(r') (r - r')^ (2) 

where d is the dimension of the embedding space. The initial probability distribution is 
assumed to be uniform, -Po(r') = KT^, with V^ the volume of the domain. 

Of particular interest is the short-time behaviour of the MSD. Let us define r the 
relaxation time of the system. Simple dimensional analysis shows that r scales as the 
square of the typical size of the domain divided by the diffusion coefficient. It appears 
intuitively clear that for t <^ r, the particle has not enough time to feel the influence 
of the boundary ^T]. The expected behaviour for the MSD is then ((5r^(t)) ~ 2dDt, 
but explicit derivation of this limit is however not always straightforward. Indeed, the 
solution of the diffusion equation is usually written as an infinite series, whose coefficients 
are related to the eigenvalues of the Laplacian. The latter actually depend on space 
dimensionality and on the particular shape of the domain. Unfortunately, an analytical 
expression for the eigenvalues is in most circumstances not available — with the notable 
exception of the problem in ID. In this article, we focus on 2D diffusion in a circular 
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domain and 3D diffusion inside a sphere. In both cases, the MSD is written as a sum of 
reciprocal powers of zeros of derivatives of Bessel functions. The goal of this note is to 
derive explicitly the short-time limit of the MSD. Though this point might sound rather 
academic — since the result is known a priori — , it certainly deserves some comments. 
To the best of the author's knowledge, this questions has never been raised before. The 
derivation outlined in this article is based on the expansion of entire functions as an 
infinite product [12]. This approach, originally developed by Euler and Rayleigh jTHj, 
has been extended recently by Muldoon and collaborators to derive convolution formulas 
for Rayleigh functions ^J E] . Adapting this method to the relevant Bessel (or related) 
functions allows us to recover the short-time limit of the MSD. The remaining of the 
paper is organized as follows. We first recall in Section|21some classical results concerning 
ID diffusion. Then we present the calculation of the MSD on a disk in section 121 and 
inside a sphere in section |31 Some concluding remarks are finally drawn in the last 
section. In order to be complete, we give in the Appendix the corresponding Laplace 
transforms of the MSD. 

2. Diffusion on a line 

To illustrate our point, we first consider ID diffusion restricted to a segment of length 
L = 2a. The Green's function G {x, x', t) satisfies the diffusion equation (^ together 
with the boundary conditions 

^(0,^',t) = ^(4^',t)=0. (3) 

Defining tlie relaxation time r = a^/D^ the solution is written as a series jlbj 

1 2 V — ^ / TlTt \ ^ t / T17TX \ ( TITCX \ 

G{x,x',t) = - + -2_^exp -[-YJ r ^°H~^J^°A~^/ ' ^^"^ 

71=1 '- -' ^ ' 

Accordingly, the MSD {6x^{t)) = {{x{t) - x(0))2) reads 

-L j-L 



{Sx\t)) = - dx dx'G{x,x',t){x-xy 




oo r 2 



96 v:^ 

p=0 



- {2p+lf- 
4 r 



(2p+l) 



(5) 



As expected, it saturates to {6x'^{t)) = L^/6 in the long-time limit t ^ t. More 
interesting is the behaviour of the MSD at short times. Indeed, writing the Taylor 
expansion of the exponential up to first order, we get 

, L^ 16L'^ 1 16Dt ^ 1 ft^\ 

^^^^'^^- 6 n- Z.(2p + 1)4+ vr^ L. ^2p + ly ^ ^ [r^ J ' 

From the well-known results of the series 
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we find that the MSD does not depend on the size L of the domain at short times, and 
recover in this hmit the celebrated result 

{Sx\t)) =2Dt t<^T . (6) 



3. Confined diffusion in a circular domain 



3.1. Green's function in 2D 

If the previous series are easily evaluated, the corresponding problem in two dimensions 
is more involved. In this section, we focus on the random motion of a tracer particle in 
a 2D circular domain of radius a. In this geometry, the Laplacian reads 

2 _ 1^ _d_ ^^ 
p dp dp p^ d^p^ 
where r = (p, Lp) are the coordinates of the particle. The eigenfunctions of the radial part 
of V^ are the Bessel functions of the first kind J„, which are solutions of the following 
differential equation [12] 

xV;'(x) + xJ'^{x) + (x^ - n^)Jn{x) = . (8) 

The reflecting boundary conditions on the frontiers of the domain involve vanishing 
radial derivatives. It can be shown that for n G N, J/^ has an infinite number of real 
zeros, all of which are simple with the possible exception of x = P^ . In what follows, 
we shall note a„m > the m*'' positive root, J^(a„m) = 0, the zeros being arranged 
in ascending order of magnitude: < an\ < an2 < . . .. The Green's function of the 
problem is then given by the following expression [TBj 
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where we define the time-scale r 



a? /D. From this solution, we obtain the MSD 
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(10) 



with 5r{t) = r(t) — r(0). As expected, it saturates to {6r'^{t)) = o? in the long-time limit 
t^T. However, the short-time limit t <^ r is more difficult to obtain. Writing the first 
terms of the Taylor expansion, we get 



{ir\t)) 



m=l 



"L 
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Im 
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wtY^ 



m=l 
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o r^ 



Obviously, the short-time behaviour should depend neither on the size nor on the 
particular shape of the domain. Rather, one expects to find {6r'^{t)) ~ 4Dt for t <^ t. 
We therefore have to get the numerical value of those sums in order to definitely check 
this point, and we now focus on generalized Rayleigh functions in order to elucidate this 
question. 
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3.2. Generalized Rayleigh functions 

At first sight, the evaluation of the infinite series seems rather difficult to achieve. But 
as we shall see here, an explicit expression for the zeros aik is actually not required. 
We remind that the ai^'s are zeros of the Bessel function Ji{x), the latter being given 

bynnmi! 

n=0 ^ ' k=l ^ l*^^ 

The validity of the infinite product expansion follows from general properties of entire 
functions [12]. Let us then define the generalized Rayleigh function | 



5n = 5^(au.)"'" (12) 

fc=i 

with 72 > 1 an integer. Notice that the series is always well defined since aik ~ (A; — jjn 

for large values of k ^\. To evaluate 5*^, we follow a method originally developed 

by Euler and Raleigh ^21 ^1 El- At a first step, we introduce a new function 

f{x) = J[(^/x), with X > 0. The zeros of / are C,k = alf.. From the representation 

of J[ as an infinite product (fTTj) . we write / as 



^<^'^^n(-i) 



fc=i 

Next, we may differentiate this expression logarithmically to obtain 
fix) ^ 1/6 



^ 1 - x/a ^ £^ ^ V a 



/(^) tt^-^/^^ tt^^^va. 

The last series is absolutely converging for < x < 1 (since aik > 1 fTT"). As a 
consequence, the order of summation can be interchanged. Also evaluating the left- 
hand side of the equality, we finally get 



n=l X \ V / 

All the required quantities can be deduced from this result. First, we evaluate the 
series ^^j^ Sn- This is done easily if we choose x = 1 in (jTHl), and together with (jHl) we 
find 

ra=l 

Next, each series Sn can be obtained by identification of the coefficient of x". Given the 
Taylor expansion (fTTj) of the Bessel function Ji, we readily get 

^1 = ^ • (15) 

Though the series with n >2 can be evaluated following the same procedure, they will 
not be required for our purpose. 

I Original Rayleigh functions are defined as the sum of reciprocal powers of zeros of Bessel functions. 
Here, we consider zeros of derivatives of Bessel functions. 
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3. 3. Statistical properties 

Coming back to the original problem, we first need to calculate the series 



oo oo 

2n 



2^ («2 _i) = i^ ^ i-a-2 =2^2^ ("1™ 

„^1 V"lm -^; ^^1 "Im -L "Im m=l n=l 

Interchanging the order of summation, we obtain 



E7jx-n) = Es- = 5- (16) 

m=l ^"l™ ^^ n=l ^ 



The other term is not more difficult to evaluate. Indeed, we have 



oo 



m=l J^'" ^ J^'" ' m=l ^ ^'" ' m=l ^"^ ra=l 

SO that we find 

°° 1 1 

£l «L («L - 1) ^ 8 ■ ^^^^ 

We indeed recover the right behaviour for the MSD at short times 
{5r\t)) = a^-8a^x^ + 8Dtx^ + o(^ 

= ADt t<^T (18) 

as expected for 2D diffusion. 

At this point, it should be noticed that the method cannot be extended to 
systematically get each order of the expansion of the MSD. Indeed, since the 
corresponding series do not converge, it is not possible to interchange the order of 
summation with the Taylor expansion. Actually, this remark already holds for the 
problem in ID — see ©, and reflects the fact the MSD is a non- analytical function of 
time. Indeed, it can be checked from the large s expansion of the Laplace transform ()A.2|1 
of the MSD that the next to leading term is proportional to t^/^. 

4. Diffusion confined inside a sphere 

4.I. Green's function in 3D and statistical properties 

Similar questions arise for the diffusion of a particle confined in a sphere of radius a. In 
this geometry, the Laplacian reads 

^2 1 ^ 2 "9 1 d . ^ 1 d^ 

V' = -Tr^V + ^~^7iH^^^^+ 2 ■ 2n 7^- (19) 

Eigenf unctions of the angular part of V^ are the spherical harmonics Y;m(6', </>), with 
eigenvalues — /(/ + 1). Eigenfunctions of the radial part are the spherical Bessel functions 
defined as ji{x) = {tc/2Y/'^x~^^'^Ji+i/2{x), with / a positive integer. They satisfy the 
following differential equation [T^] 

x''j'/ix) + 2x][{x) + {x^ - l{l - l))3i{x) = . (20) 
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We also need to enforce reflecting boundary conditions on the surface of the sphere. 
To this aim, we define /3;„ the (non-zero) zeros of the derivatives of the spherical 
Bessel functions, jKPin) = 0. They are arranged in ascending order of magnitude: 
< Pii < f3i2 < ■ ■ ■■ The normalization condition for the spherical Bessel functions 
is [E] 



drr^ji (pin-) 
^ a/ 



2AI 



(/3f„ -/(/ + !)) J, (An)' 



so that the Green's function for diffusion confined inside a sphere reads 

1=0 m=-l 

fl jiyPln-jjiyPln-^ 



X 



oo n2 

Pin 

/32 - 



n=l 



^L 
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exp 
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Here, we used the same definition r = o? / D. 



(21) 



(22) 



^.£ Statistical properties 

From (I22I), we readily find (r2(t)) = 3a75. The MSD is then 

(^r^t)) = — -2(r(t)r(0)cos7) 



with 7 the angle between the vectors r(t) and r(0). Rewriting C0S7 = cos6'cos6'' + 
sin 9 sin 9' cos(0 — 0') in terms of spherical harmonics 

1 



47r 



cos7=^ E n*™(^,0)nJ^',0') 



m=— 1 



one can easily convince oneself that only the terms with / 
contribution to the mean value. We get after some algebra 



(23) 
1 give a non-zero 



[5r\t)) 



6a2 



t 



In 



T 



1 



PL iP' 
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00 

12a^^exp -0-_ 

n=l 

As expected, the MSD saturates to {6r^{t)) oc a^ in the long-time limit t ':$> t, whereas 
for short times t ^ r we find 

^2 



fi^2 ~ 1 



^ /32„ (/32„ - 2) 



12DtY^ 



ra=l 



/5fn - 2 



-'^ 



The situation is similar to that in 2 dimensions, and we now detail the calculation of 
the infinite series. 
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4-3. Short-time behaviour 

First, we rewrite the series that appear in the Taylor expansion of the MSD. After some 
elementary manipulations, we find 



CXD oo / /TT \ ^" _ oo 



2^ /92. _ 9 9 2^ 2^ /9, , 9 2^ 



— /?2 2 ~ 2 ^ ^ I /?ifc / 2 — 

fc=l '^1'= fe=l n=l V^'i^/ „=l 

where the series o"„, with n an integer > 1, is defined as 

V 



<Jn 



k=l 



Note that /^n > \/2 17 so that all sums do converge. The other term is given by 

Ei i ^^ — \ i i ^^ — \ i 



k 

oo 




The central point is to evaluate the series cr„. To this aim, it appears convinient to 
consider the function j[(v^2x), whose zeros are given by Cfc = Pikf^- W^e then write 
j[{^/2x) as an infinite product § 

oo . 



where the normalization factor comes from j[(0) = |. Taking the logarithmic derivative, 
we are directly lead to 



E-'^-" = -a/!-^^- (26) 



From this result, we may draw the following conclusions: 

• for X = 1, the numerical value of the right-hand side can be obtained from 
the differential equation ^K^ satisfied by the spherical Bessel function. Since 
2j'{{V2) + 2V2j[{V2) = 0, it is found 



oo 



Ea„ = l. (27) 



n=l 



the series o"i is obtained by identification of the linear term in (j26|) . The Taylor 

30 



expansion of the modified Bessel function being ji{x) = I — In + 0(^^)) we obtain 



ai = ^ . (28) 

Bringing everything together, we obtain the expected result for 3D diffusion, namely 

{6At)) = ^ - 12a^ X 1 (l - + 12Dt x ^ + O (^ 

= 6Dt t<^T . (29) 

§ Compare with expression (11) of ^1] with a = 0, 6 = 2, c = —1, and i^ = 3/2. 
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5. Conclusion 

To summarize, we have derived explicitly the mean square displacement for confined 
diffusion in a disk or a sphere. Since the results involve zeros of derivatives of Bessel 
functions, the numerical values of the relevant series are not easily obtained. To by-pass 
this difficulty, we used a powerful method based on the expansion of entire functions 
as an infinite product. This allowed us to explicitly check that {6r'^(t)) ~ 2dDt in 
the short-time limit. Extension to more general domain shapes or space dimensions 
is straightforward, provided that the eigenfunctions of the Laplacian are known. This 
approach can also be adapted for adsorbing or mixed boundary condition, and might 
as well be relevant for the study of polymer chains under confinement. Indeed, the 
partition function of a polymer satisfies a diffusion- like equation jTHj, where the time 
variable is replaced by the polymerization index N. In this case, the limit t <^t would 
correspond to Rg = ^/nWJQ <ti a, with a the typical dimension of the cavity and h the 
Kuhn length. 

Finally, we have to mention that the actual problem becomes more involved in 
the (highly relevant) physical situation of a confined colloidal suspension. Indeed, 
hydrodynamic interactions with the bounding surface have to be accounted for if the 
particles are moving in a liquid, and the diffusion coefficient is not expected to be 
constant anymore. Instead, theoretical jTH] and experimental [201 1^ I22j works have 
revealed that the mobility is in fact a function of position — or, more precisely, of the 
distance to the surface. In addition, the diffusion coefficient may also be time- dependent 
when the bounding surface is soft, like a liquid-liquid interface or a fiuid membrane. In 
some situations, it may even be possible to observe anomalous diffusion, in which case 
the mean square displacement follows a power law ((5r^(t)) ~ t", with a 7^ 1, over a 
notable range of displacements and observation times [231 121] • 
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Appendix A. Laplace transform of the MSD 

In this Appendix, we give the expression of the Laplace transform of the MSD. 
Defining ad{t) = {Sr^{t)), where d is the dimensionality of space, the Laplace transform 
ad{s) = /q°° dtexp[— st]crrf(t) is evaluated by solving directly the diffusion equation in 
terms of the variable s rather than t. Average values are then easily evaluated since no 
infinite series are involved. Notice however that ad{s) can also be obtained by direct 
Laplace transform of the results ©, (fTTH) and (j^ . though the method described in this 
article has to be extended in order to evaluate the corresponding series. The result is 



^ 2D 



S2 



tanh (y/sr) 



(A.l) 
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for diffusion confined on a segment. For diffusion in a disk, one gets 
^ , , AD 



^_ f hiV^) 



with /i the modified Bessel function. Finally, the Laplace transform of the MSD for 
diffusion confined inside a sphere is 



6D 



^3 S = ^- 



tanh (i/ir) ~ V^t 



(A.3) 



{sT + 2) tanh (a/st^) - "^y/sr 
Asymptotic limits follow in a straightforward way. 
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